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ABSTRACT 

This is the written version of a set of lectures on perturbative QCD that were de- 
hvered to a mixed audience of young theorists and experimentahsts in the course of the 
XXII International Meeting on Fundamental Physics. These notes are virtually a verba- 
tim transcription of the lectures. The selection of topics is somewhat arbitrary, but two 
basic points are emphasized: the rationale behind QCD and how ongoing experiments, 
such as those taking place in LEP and HERA, contribute to our understanding of strong 
interactions. The lectures are thus roughly divided into three parts: foundations of QCD, 
determination of as and jet physics, and deep inelastic experiments. 



^ Lectures delivered at the XXII International Meeting on Fundamental Physics "The 
Standard Model and Beyond" , Jaca, February 1994. To be published in the proceedings. 
* E-mail: espriu@ebubecml.bitnet, espriu@greta.ecm.ub.es 
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1. - Introduction 

It is not the purpose of these notes to give a detailed account of Quantum Chromodynamics, 
the theory of strong interactions. There are many good textbooks^-* that are exceUent at 
doing this job. We do not pretend to treat exhaustively any particular aspect of Quantum 
Chromodynamics or attempt to provide a general review of the status of this theory, either. 

These lectures are essentially addressed to young experimentalists without too much 
theoretical background in Quantum Field Theory. We have attempted to single out 
amongst the theoretical foundations of QCD the ones where Quantum Chromodynam- 
ics really hinges on, trying to present them in terms as simple and physical as possible. 
Then we move to the connection between theory and experiment. This is not a simple 
problem in QCD. It is well known that the fields and particles we know how to compute 
with (with the simplest tool at our disposal, perturbation theory) are not those that are 
observed by experimentalists in their detectors due to the phenomenon of confinement. 
Quarks and gluons are real, but they cannot be detected as free particles as they are be- 
lieved to be confined inside hadrons. In view of this is quite remarkable that there are 
theoretical techniques enabling us to put the theory to very stringent tests. To mention 
just a significative measurement, the value Q!s(M|) = 0.123±0.006 from an analysis of the 
process e+e" — > qq has been reported in this meeting^). This is an impressive accuracy 
for a number that is not, strictly speaking, observable due to confinement. Even more so 
when we compare this precision to the one we had only a few years ago. Reaching this 
level of accuracy is one of the achievements of LEP. 

Perturbative QCD can be applied to inclusive processes provided that the character- 
istic momentum transfer is large enough. In these lectures we will discuss is some detail 
these techniques and apply them to the determination of ctg through Rhad Rt, or from 
jet topology. Another test of perturbative QCD is, of course, deep inelastic scattering. 
Here the commissioning of HERA has opened a new kinematical region where it will be 
possible to study the onset of non-perturbative effects. The exploration of this region is a 
fascinating subject interesting on its own right. 

Due to the lack of space we have not included two sections that, in our view, should 
be in any general review of QCD. The first one concerns the so-called "spin of the proton" 
problem^). Another topic that is not covered at all is the study of the photon structure 
functions. We direct the interested reader to ref. . 

We have not really made any serious attempt to provide a complete set of references. 
Those provided merely reflect personal tastes. A more complete bibliography can be found 
in some of the references, as indicated. 

2. - Why QCD 

Since Dirac we know that the vacuum in Quantum Field Theory is a rather peculiar object. 
The vacuum is defined by the property that all its negative energy states are filled, while 
those with positive energy are empty. A virtual photon with sufficient energy will promote 
a negative state to a positive one, thus creating a particle-hole (electron-positron) pair. 
Eventually, the positive energy state will decay filling again the hole in the Dirac sea with 
the emission of a off'-shell photon which, in turn, will excite another pair and so on. This 
is illustrated in fig. 1. 
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Fig. 1.- The Dirac sea. 



In the quantum world this process takes place virtually all the time. The uncertainty- 
relation AEAt > h allows that for a short time a pair is created, then annihilated. The 
vacuum is thus constantly populated by virtual electron-positron pairs, becoming a dielec- 
tric medium at quantum scales. A medium that partially screens charges. The screening 
is more and more effective at larger distances since more and more virtual pairs enter the 
game. Therefore 

« ^ «e//(r) 

ro is some scale at which we choose to measure a. At large scales the effective electro- 
magnetic coupling decreases, becoming larger at short distances. While it requires some 
work to get the detailed form of eq. 1, the reasoning behind the behaviour of the effective 
charge for r ^ tq is so general that it is difficult to see how it could go otherwise. 

Yet, the opposite behaviour takes place in Quantum Chromodynamics. Fig. 2 shows 
the elastic and total cross sections for tt scattering off protons. The cross sections show 
a rather complicate resonant behaviour at low energies but become much simpler at high 
energies. This suggests that strong interactions become stronger at low energies, i.e. long 
distances. 



Fig. 2.- A typical hadronic cross-section at low energies^-* 

In a completely different regime deep inelastic scattering (fig. 3) showed 25 years ago 
the existence of hard processes inside nucleons, at very short distances. One can get a 
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rough picture of these processes by just assuming that quarks are approximately free if 
one looks at distances <^ 1 (GeV)~^. 



Fig. 3.- The kinematics of Deep Inelastic Scattering. 

One needs a rather peculiar type of theory. It must be a theory whose elementary fields 
are fermions (quarks), but weakly interacting at short distances. In addition we would like 
it to be a renormalizable theory to be able to apply the full machinery Quantum Field 
Theory and this pretty much forces us to work with gauge theories^-* . Trivial modifications 
of QED will not work because we need anti-screenning of charges rather than screening. 

Several hints as to the way to go come from the successful quark model. We have 
for instance the A"'""'" = \u^u^u^). Being the lightest hadron with this quark contents we 
expect to have the three quarks in the ground state, hence in a symmetric wave function. 
This is in contradiction with Fermi statistics. The contradiction can be solved if we admit 
the existence of a new quantum number a and 

I^T^T^T) = J-e<-^'^\ululul). (2) 



Indications in the same direction also come from two very different processes. One is the 
cross-section for e+e" — > hadrons. Assuming a free quark model (which as we mentioned 
should be good as the momentum transfer is large) we have 

a{e+e- ^ hadrons) /Q^ 



i=l 



The index i runs over all degrees of freedom coupling to the intermediate photon (with 
charge Qi). 

n = 3 u,d,s R = - 

' ' 3 

n = 4 u,d,s,c R=— (4) 

5 

n = 6 u,d,s,c,b,t R=-. 

Experimentally all these values are off by about a factor 3. The quark model can again be 
reconciled with experiment if we admit that the new quantum number can take 3 values 
and is blind to electromagnetic interactions, q'^, a = 1, 2, 3. 
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Fig. 4.- The decay 7r° — > 77. 

The other process relevant at this point is the celebrated 77 decay. This process 

goes dominantly through the triangle diagram shown in fig. 4, with a closed loop of quarks. 
This diagram is also called 'the anomaly', for reasons that we will discuss in a while. The 
calculation gives^^ 

^= 6^ fr 

Experimentally F = 7.37 ± 0.5 eV. The result is off by a factor 9 = 3^ which, again, is 
understood if we accept that additional degrees of freedom, invisible to both photon and 
pion, exist. So we learn that ordinary hadrons seem not to carry any color. Similar analysis 
show that color does not couple to the W or Z bosons either. 

The guess is now more or less obvious. Let's introduce a gauge symmetry acting on 
the new degree of freedom just as Electromagnetism acts on the electric charge with a 
gauge group U{1). It must be a new gauge group, since color does not couple to known 
gauge fields. We want a group with irreducible representations of dimension 3, for the 
quark model to fit in. Obvious candidates are 0(3), SU{2), U{2), U{3), SU{3), SO{3) 
and Sp{2). Groups such as U{2), U{3) and 0(3) can be discarded right away because 
ga/37 j^Q^ invariant tensor (and we need that for the A"'"''). On the other hand 
Sp{2) ~ SO{3) = SU{2)/Z2- Neither of these have complex representations of dimension 
3 (hence they would lead to diquark states). We are left with SU(3). 

3.- Lagrangian and Symmetries 

The QCD lagrangian is 

^ n n 

where 

= d^.W^ - d^Wl + gfa,>cWlWt, (7) 

D^^p = 8^^d^-igTl^Wl. (8) 

The generators T" are related to the well-known Gell-Mann matrices and act on the fun- 
damental representation of SU (3) 

7^a/3 = ^ {T\T''\=ij,^,T- (9) 
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With the structure constants fabc one constructs the generators of the adjoint representa- 
tion T^"^ = ifabc- 

^QCD has a local gauge invariance. If G{x) is a SU(2>) matrix, the transformation 

^l){x) G{x)ij{x) Wf,{x) G{x)W^{x)G{x)-^ + ^df,G{x)G{x)-^ (10) 

leaves Cqcd invariant. This symmetry is crucial to remove two of the four degrees of 
freedom in the field W^, {p, = 0, 1, 2, 3). 

To quantize the theory one must select a gauge. The bilinear part in the field is 

iw^;(/cV'' - k^'k'')W^ = ^WI^Mf'W^; (11) 
M'^'^ cannot be inverted to find the propagator. The way out is to add the piece 

^id^w;r. (12) 

Then 

M^"' = k'^g^"' - (1 - -)k^'k'' {M-^Y" = ^ ILJs^^ (13) 

We can now write Feynman diagrams. The added term eq. 12 breaks the local gauge 
symmetry which, generally speaking, is only recovered for jS-matrix elements. 



Fig. 5.- Interaction vertices in QCD. Only the color factors are shown. 

^QCD has several types of interaction vertices. They are shown in fig. 5. Let us now 
consider the process qq —>■ gg. At tree level the appropriate diagrams are shown in fig. 6. 



Fig. 6.- Tree level diagrams for qq — > gg. 

(c) is absent in the analogous QED process e+e~ — > 77. Due to this diagram it turns 
out that the above process has a bad high-energy behaviour when we sum over final state 
polarizations covariantly 

P =\9^.^.9^..uJ''''Kr''^^)\ (14) 
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but it is just fine if we include transverse polarizations only 

The tensors g^i, and r^jy are obtained by summing the polarization vectors e^{a)ei,{a) 
over all polarizations or over physical ones only, respectively. If we insist in keeping a 
covariant formalism in which we sum over all four polarizations something must cancel 
this undesirable high-energy behaviour. This is accomplished by adding to Cqcd the 
piece 

-d^(paDf^, (16) 

with being the same as D'^^ in eq. 8, but with the generators of the fundamental 
representation of SU{?>) replaced by those of the adjoint one. The fields have boson-like 
couplings, but are defined to have Fermi statistics. They contribute with a (-1) factor to 
the cross-section. They are not required in abelian theories like QED, but are crucial in 
QCD. 



Fig. 7.- Ghost contribution to qq — > gg cross-section. 

If all polarizations, physical and unphysical, are summed over one must accept that ip states 
can be produced even if they are ghost states with unphysical statistics. The contribution 
from the piece we have added to the lagrangian, eq. 16, is just right to reproduce the same 
results we would get keeping the physical polarizations only. In practice, in internal loops 
we have really no choice but to keep the covariant sum over polarizations and ghosts have 
to be included there. It is also possible to derive the need for the introduction of ghosts 
from more formal arguments^^, but that would take us too far afield. 

In addition to the local gauge symmetry we also have exact or approximate global 
symmetries in Cqcd- The lagrangian is invariant under the global transformation 

il){x) ex.p{-i9I)iij{x), (17) 

leading to baryon number conservation 

B = J d^xJo. (18) 

If all quark masses are equal there are additional symmetries and conserved currents^) 

ijix) exp(-irT")V'(a;); (19) 
= 7/i7^T>. (20) 
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i/j now represents a column vector containing all flavours and is an SU{Nf) generator. 
Furthermore, if all quark masses vanish Cqcd is also invariant under 



V^(x) ^ exp(-irT"75)V(a;). (21) 

= V^TmTsT^ (22) 

are the corresponding conserved currents. In the real world quark masses are not equal, 
let alone zero. The two latter symmetries eqs. 19, 21 are only approximate and this only 
for light quarks. Therefore the hadronic world is approximately invariant under U{1) x 
SU{3)v X SU{3)a- U{1) is always exact and SU{3)v is nothing but the vintage SU{3) of 
Gell-Mann-'^^) , which led to the quark model thirty years ago. 
J^QCD is also invariant under 

i;{x) ^ exp(-i^/75)V(a;) (23) 

The 'conserved' current is 

4 = V^TmtV- (24) 

However, when one computes Green functions with insertions of the divergence of the 
above current, d^J^, one gets non-zero answers. The culprit is the triangle diagram, the 
'anomaly' (the same that we met in 7r° decay, but with the two external photons replaced 
by gluons). In fact, a careful calculation shows that while the axial current (24) is conserved 
at tree level, quantum corrections spoil that conservation and, in fact, 

dMn'i^) = |^^^^^"^^;.^a/3- (25) 

The key point is that there is no way of estimating the divergent momentum integral that 
appears in the evaluation of fig. 4 without breaking the U(1)a symmetry, and this is not 
a point of mathematical finesse; it has far reaching consequences. The r.h.s. of eq. 25 is 
itself a total divergence d^K^. Then the charge verifies 

Q5 = y d^xdo4 = J d^xd^K" - J (fxdiJi (26) 

In perturbation theory all fields are small perturbations from the vacuum; they decay fast 
enough to infinity to be able to neglect all boundary terms in the integrals. By Gauss 
theorem the second integral on the r.h.s. of eq. 26 drops and 

Q^{t = +00) - Q^it = -00) = J d^xdf.Ki' = (27) 

In a non-abelian gauge theory such as QCD there are, however, some non-perturbative 
gauge configurations^^'' (i.e. configurations which are not a superposition of states with a 
finite number of quarks and gluons) that do not possess the nice long distance behaviour 
that is required for eq. 27 to vanish. These make Q5 7^ 0. A conserved charge can still be 
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defined by the integral of Jq — Kq, but it is not gauge invariant (the divergence of is 
invariant, but not i^T^ itself). There is no way to have a conserved, gauge invariant axial 
charge in QCD. U{1)a is not a symmetry of the theory. 

At this point we should probably make some contact with ongoing experiments. Can 
we test in some simple way that the color structure due to the non abelian theory is indeed 
correct? We have to take into account fig. 5 and recall some simple group relations and 
definitions 

J2 = CpSaj, Yl ^bcT^d = CASbd, TriT-T'^) = TSab- (28) 

0/3 ac 

In QCD Cf = 4/3, Ca = 3 and for generators in the fundamental representation T = 
Tp = 1/2, while in the adjoint representation T = Ta = 3. The factors Cf,Ca,Tf and 
Ta appear in jet counting rules, as shown in fig. 8 



Fig. 8.- Lowest order jet cross sections. 

As evidenced by fig. 9 the experimental agreement between the LEP data and the QCD 
predictions is perfect. 



Fig. 9.- QCD color factors as measured by ALEPH and DELPHI. From^^) 
4.- Renormalization 

Beyond tree level most Feynman diagrams are ultraviolet divergent. Take for instance the 
first diagram of those contributing at one loop to the gluon propagator (fig. 10) 
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Fig. 10.- One-loop contributions to the gluon self-energy. 

Neglecting external momenta, the integral over the momenta of the internal particles is of 
the form 

To make sense of the theory and get a finite result we must introduce a cut-off A and 
counterterms. A typical method is to perform a subtraction at some (f' = —jJ,^. For 
instance, for the self-energy of the gluon propagator 

U{q^) - U{-fi^) = Uuiq'') = finite. (30) 



Fig. 11.- Adding counterterms to make Green functions finite. 

Alternatively we can make sense of the integrals using dimensional regularization by con- 
tinuing the dimensionality from 4 to n = 4 -|- 2e, 

/d'^k f d^k 

(2^^7(2^' ^^^^ 

and subtract just the poles in 1/e (minimal subtraction, MS) or also the — log47r 
that always accompanies the singularity in 1/e (improved minimal subtraction, MS). For 
instance, for the second diagram in fig. 10, namely the quark contribution to the gluon self- 
energy, one has the following result after computing the integral in n = 4 -|- 2e dimensions 

nfa^) = -^Mj + 7. + log ^ + ...). (32) 



Using the MS and MS schemes one gets 

nMs(g') = -^5a,{iE + log ^ + • • •) n^(5') = -£^«''(i°g ^ + • • ^^^^ 



Admittedly this looks totally adhoc. It would seem that we can get rid of any divergent 
integral by just adding the appropriate counterterms. Of course, since these are sort 
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of arbitrary, we would then be able to obtain any result we want. This is not so, of 

course. A renormalizable theory is one in which the necessary countcrtcrms to all orders 
in perturbation theory are generated by redefining the fields and paranietres of the original 
lagrangian and nothing else, and this is highly non-trivial. QCD is such a theory: if we 
redefine 

_ 7 -7-3/2 _ ^ ^-1^-1/2 _ ry 7-1/2^-1 _ ^1/2^-1 

9 — ^iym^sym9 — ^1^3 ^zym9 — ^if^zym^2f9 — ^5 ^zym9i 

= ~Zy\. (34) 

it is possible by a suitable non-unique choice of the Z's to make QCD finite. Since all Z's 
are of the form Z = \ -V with AZ beginning at 0{g'^) one must replace the original 
^QCD one started with by Cqcd + ^^qcd- ^^^qcd contains the necessary counterterms 
that must be added. Equivalently, we can work directly with the QCD lagangian with all 
fields and parametres replaced by the 'bare' ones, defined through eq. 34. Observables 
will be obtained by the prescription (A being the ultraviolet cut-off) 

r(Pi,i/,/i) = lim Z~yl^Z~p^'^VQ{pi,gQ,K). (35) 

A— >oo 

On the l.h.s. we have renormalized Green functions or amplitudes expressed as a function 
of renormalized parametres. On the r.h.s we have bare Green functions or amplitudes 
expressed as a function of bare parametres. In dimensional regularization 

r(Pi, g, IJ) = lim Z~^jl^Z~p^'^VQ{pi, go, e). (36) 

e— >0 

fi is the subtraction scale and m and k are the number of gluon and quark external lines, 
respectively. The effect of counterterms is to replace the dependence on the cut-off (A, e, ...) 
by a dependence on /i. Unlike in QED where the natural scale is mp, or the Electroweak 
theory where the natural scale is M^, there is really no preferred way of choosing the 
counterterms. The only requirement is the fulfillment of the Ward identities'^-* . In practice 
MS and MS are the most useful, particularly the latter that seems to lead to perturbative 
series with a faster convergence. In the MS scheme all renormalization constants are just 
poles in e 

^ = 1 + — - + { — f{^ + -+... = 1 + ^ _ _ _ 37 



5.- Renormalization-group Equations 

Note that regulated quantities depend on a cut-off (A, e, ...) and that the renormalization 
of fields and constants through eq. 34 trades the dependence on the cut-off by some scale 
[I. Yet, physics cannot depend on at all. If you change you must change at the same 
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time the value of your renormalized parametres to make up for the difference. A simple 
way to encode this observation is the following. Let's write 

T{pi, as, ij) = Zr(/x, e)ro(pi, q;°, e) (38) 

To is obviously independent of the subtraction scale. Therefore 

(We neglect everywhere the dependence on the gauge parameter, as well as the quark 
masses. Of course they have to be properly taken into account. Physical on-shell ampli- 
tudes are ^ independent.) From eq. 39 

or, defining the so-called /3-function and the Green function anomalous dimension 7r, 

l^-^ = asP{as) y"^^r = 7r, (41) 

d d 

{li— + (3as^ 7r)r(pi,as,//) = 0. (42) 

oil oag 

This is called the renormalization-group equation -^^^ . Let's investigate its consequences. 
We can always write F as a function of dimensionless variables by pulling out ji^ {D: 
dimensionality of F) 

T{Xpi,as,fJ,) = ^^F(^^^,a,). (43) 

Hence 

(A^ + /X— - D)r{Xp,, as, fi) = 0. (44) 
Using the renormalization-group equation we get 

{-- +I3as^ 7r + D)r{e'p„ as, /i) = 0. (45) 

ot oag 

Prom subtraction scale independence arguments we have been able to establish an equation 
concerning the dependence on the external momenta. We can formally solve this equation 

F(eVi, asifJ'), fJ') = ex.p[tD - [ dt-fr{as{t))] x T{pi, as{t), fJ,), (46) 

Jo 

as{t) is just as{e*iJ,), i.e. the same coupling constant but renormalized at a different scale. 
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From eq. 46 we see that when we scale the external momenta in a Green function or 
amplitude the change is absorbed 

(i) in a multiplicative factor that depends on the anomalous dimension as well as the 
engineering dimension of the amplitude and 

(ii) in a redefinition of the coupling 

asifi) as{e*ii). (47) 



The renormalization- group evolution of the parametres in the theory (in this case exem- 
plified by the coupling constant ctg) governs the scaling behaviour. We have to find which 
is the evolution of as under a change of /x. This is actually a very simple question. We 
just have to solve the equation 

= c^s/3{as). (48) 

Since the n dependence of as is introduced via counterterms, To compute P we have just 
to find the relevant renormalization constants (see eq. 34). 

Of course (3 is evaluated in perturbation theory and, accordingly, the differential equa- 
tion eq. 48 is also solved in perturbation theory. The solution will only make sense as long 
as the expansion parametre is small. At one loop the solution of eq. 48 is 

c^.(eV)= , ,y (49) 

1 - -/3ias{i^)t 

In QCD Pi (the first coefficient of the function) is 

(recall that in QED Pi — 2/3). Note that Pi is negative if Nf < 16. If Pi is negative, at 
larger momentum transfers, where the relevant scale will be e^/j, and not /i, as will actually 
decrease. The solution of the renormalization-group equation will actually become better 
and better at higher energies. 

6.- Asymptotic Freedom 

Whenever we speak of 'higher' energies we are implicitly assuming the existence of some 
characteristic QCD scale. Looking at eq. 48 we note that 

t=\\og^?= I = i:ias) + a (51) 

2 J asP[as) 

C is an integration constant. Therefore t — ip{as) is a constant of motion along the 
renormalization-group trajectory. At one loop, we plug Pi in eq. 51 and get 

1 + W^TJa ^^=\ =^ «.(/^) = p,, r- (52) 

2 Piasifi) 2 f\og{fi^/AlcD) 
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If we renormalize at scales much larger than Aqcd the renormalized coupling constant will 
be small and working at one loop will be justified. Table 1 shows a number of values for 
as extracted from different LEP observables. The relevant scale here is M^, much larger 
than any hadronic scale. At such energies perturbation theory is clearly meaningful. 



Table 1.- Determinations of ag{Mz). Prom-'^^^ 
Let us recall the physical contents of eq. 46. We can write it symbolically as 



Fig. 12.- Asymptotic freedom visualized. 



Ignoring for a second the overall factor, when we scale up the momenta we have exactly the 
same amplitude, but renormalized at scale e*//, i.e. replacing q;s(//) by Q;s(e*//). Obviously 
the scaled up amplitude will correspond to a theory that interacts more weakly. In the 
t — > oo limit we will have a free theory. This is called asymptotic freedom and it is one of 
the characteristic signals of strong interactions. QCD has it, while QED has not. 

From fig. 12 we learn something else. Let us imagine that for some reason radiative 
corrections are small for some amplitude with momenta pi and the coupling constant 
renormalized at scale fi. (For instance, in QED we may decide to choose the counterterms 
and define the renormalized coupling constant in such a way that the classical formula 
for Thompson scattering is strictly valid at all orders in perturbation theory for on-shell 
particles.) The renormalization group tells us is that if we scale the external momenta and 
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the renormalization scale in the way prescribed by eq. 52, radiative corrections will also be 
small for the scaled amplitude. Each physical process has a characteristic scale. Choosing 
this scale as subtraction point typically optimizes the perturbative series. 

We expect the coupling constant ctg to be small at LEP energies because Mz 3> Aqcd, 
but which evidence do we have that it actually runs according to eq. 49? Well, actually 
as is not small enough for two-loop effects to be completely neglected so to answer this 
question in a precise way we have to work a bit harder and compute the two-loop /3- 
f unction. The result is 

/52 = -| + gA^/. (53) 

Incidentally, both Pi and P2 are scheme independent (but not /?3 and beyond). Note that, 
for Nf < 8, /32 is also negative. Repeating the steps that led to eq. 52 we get at the 
two-loop level 



1 



log/x^ + 



1 



+ -2 logas(//) + 0{aa) = -logA 



2 

QCD: 



(54) 



whose solution is 



127r 



(33 - 2Nf) log(^VA^c^) 



1-3 



153 - IQNf loglog(/xVA|gj^) 
(33-2iV^)2 1 log(/.VA|^^) 



(55) 



Working in the MS or MS scheme changes the value of a.s{p) but the difference is 0{a1). 
Therefore the numerical value of Aqcd does actually depend on the renormalization scheme 
one is using (one does not see this at the one-loop level precision). For instance, at two 
loops Ams and Aj^ are related through 



A2 



A?—. 

MS 



(56) 



Eq. 55 is thus our starting point to check whether the behaviour predicted by the renor- 
malization group has factual support. However, there is one more thing that we have to 
account for when we choose to work in the MS or MS schemes. Decoupling of heavy 
fermions is not manifest in any of these schemes. In practice one works with the number 
of quarks that are excited at the energy one is. For instance, at scales well below and well 
above the charm threshold we have, respectively (at the one-loop level) 



12- 



271og(/xVA2(3)) 



asifJ,) 



12- 



251og(//VA2(4))" 



(57) 



The coupling constant is obviously continuous as we cross the threshold, but Aqcd is not. 
Demanding continuity on as leads to the matching condition 



A(4) 



A(3)\^ 



rrir 



A(3). 



(58) 
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At two-loops this gets modified to-^^^ 



If one insists in keeping A(3) beyond the charm threshold the corrections will be large; 
perturbation theory breaks down. 

The two loop evolution of the QCD coupling constant is shown in fig. 13 for different 
values of A-^^^. On top, the values obtained from different experiments at different scales 
are shown. Asymptotic freedom is very nicely confirmed. We cannot go to energies much 
lower than those shown because the coupling constant grows very quickly and thrcc-loop 
effects and beyond are relevant; in fact, perturbation theory becomes meaningless. In the 
coming sections we will discuss some of the ways of determining ctg that are shown in this 
plot. 



Fig. 13.- Evidence for next-to-leading scaling. From^^^ 
7.- Confinement 

Aqcd sets a natural scale in the theory. Well above Aqcd perturbation theory makes 
sense. Of course perturbative QCD at large enough energies describes a world of quasi-free 
quarks, interacting with Coulomb-like forces. We know very well that hadronic physics is a 
very different world with quarks are confined into colorless hadrons. As soon as ~ ^qcd 
perturbation theory is unreliable. It simply cannot explain confinement. 

What does confinement actually mean? One popular interpretation is that 'it is not 
possible to detect an isolated quark or gluon'. The problem with this definition of con- 
finement is that in electromagnetism, which is certainly not confining, it is not possible 
to detect an isolated electron either. Electromagnetism (as well as Quantum Chromody- 
namics) is long-range (mediated by a massless particle) and plagued (just as QCD) with 
infrared divergences. Both QED and QCD observables have to be inclusive enough. There 
is one difference, of course, and that is that photons have no U{1) charge, while gluons 
do carry SU{3) charges. Far away from the interaction point you can hope to be able 
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to measure the electric charge carried by one particle (or rather what the experimental 
resolution defines as a particle which is actually the electron surrounded by a cloud of 
soft photons), but even if you were able to construct a detector that measured color you 
probably would not be able to identify in any way the color of the quark itself. 

In fact, there is another definition of confinement that tells to you that the chances 
of actually seeing a (gluon-dressed) quark are small: 'there is a force between quarks that 
does not decrease with distance'. There is indeed phenomenological evidence (which is 
supported by lattice analysis) that the interquark potential in QCD is of the form 

V{r) ~ aAlcor - ^ + • • • (60) 

The first term is a confining quark potential. The constant a has to be ~ 1 because A.q(j^ 
is the only dimensional quantity at our disposal. The Coulombic part is called the Luscher 
term and plays a crucial role in heavy quark spectroscopy 

We like this second definition of confinement better, because the first one is far too 
imprecise. In order to see that it is useful to recall a toy example suggested by Georgi^^^. 
Imagine a world in which we tune Aqcd in such a way that Q;s(lGeV) — 1/137. Since 
the coupling constant is so small, perturbation theory works wonderfully at such energies. 
The proton would be a bound state of quarks (bound by Coulomb-like forces that is) with 
mass roughly Sm^. Its size would be dictated by the Bohr radius, about 1000 times the 
size it has in our world. The inhabitants of this world would certainly not understand the 
first definition of confinement, since the confinement radius would be ~ ~ lO^-'^cm. 
They would see quarks as you see electrons. 

Even in our world the situation is somewhat similar to that of the toy world for very 
heavy quarks. Indeed as{mt) is small (say ~ 0.1). The Bohr radius is tq ~ 10~^ fm, much 
smaller than Ag^^. The coulombic part of the interquark potential largely dominates. (At 
such short distances the linearly rising potential is not at work, the leading confinement 
effects are ~ r^, as discussed by Leutwyler some time ago-'^^^ but they can be safely 
neglected at first approximation.) 

Bottom and charm are in a somewhat intermediate position. as{mb) is still relatively 
small. The Bohr radius is 10~^ fm, smaller but comparable to Ag^^. Spectroscopy is 
basically perturbative, at least for the lowest levels, but some non-perturbative effects are 
visible. Charm is really no-man's land. Both perturbative and non-perturbative effects 
compete even for the ground state n = 1. For light quarks the Bohr radius is several fm 
and the confining potential is fully at work. 



Fig. 14.- The QCD string. 



17 



The existence of a confining potential leads to very large multiplicities and jets. The 
situation is visualized in fig. 14. One can imagine a quark-antiquark being formed at 
the primary vertex then moving apart. Part of their kinetic energy is deposited in the 
interquark potential as they move away. Very quickly a separation is reached where 
the energy deposited is enough to form a new quark-antiquark pair, 

^QCDfm - 2mq, (61) 

at that moment the quark-antiquark 'string' breaks and the process is repeated until the 
average relative momentum is small enough and hadronization takes place. 

There is a lot of physics in the string picture. We can think of color forces being 
confined in some sort of tube or string joining the two moving quarks. The chromodynamic 
energy is thus stored in a relatively small region of space-time. If this picture is correct we 
should expect hadronization to take place in this region in preference to any other. This 
is indeed the case; in three jet events (which originate from qqg, with a hard gluon) there 
is a clear enhancement of soft gluon and hadron production in the regions between color 
lines (representing the gluon by a double color line, or qq state), and a relative depletion 
in other regions. This phenomenon is called color coherence^^^ . 



8.- Rhad 

The observable 



cr(e+e hadrons) 
^had = — / I _ 1 _x \pl) 



is probably the cleanest and simplest observable in QCD. It is fully inclusive and can be 
computed through a dispersion relation which is symbollically depicted in fig. 15 



Fig 15.- The imaginary part of the self-energy describes the 7* — > hadrons width. 

The sum over intermediate states on the r.h.s. of fig. 15 runs over all hadronic states. 
We can just as well compute this sum using another resolution of the identity — the one 
that is provided to us by perturbative QCD, i.e. in terms of quarks and gluons. Rhad has 
been computed in this way up to third order in 

= fi„[l + + r,(^f + .3(^^)= + ...], (63) 



^ We shall write from now on Q!s(//) or agdi^) equivalently. We have already argued that 
the right renormalization scale should be the typical momentum transfer of the process. 
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where (f' is the typical momentum transfer, and in the MS scheme 

r2 ~ 2.0 - ^.12Nf rs = -6.637 - 1.200iV/ - 0.005iV| - 1.240^^^. (64) 

For Nf = 5, r2 ~ 1.4 and ~ -12.8. We have seen that aa{Mz) ~ 0.12. Then 

Rhad = i?o[l + 0.04 + 0.002 - 0.0008 + . . .] (65) 

The convergence of the series does not look bad, but it is not very good either. The values 

for the coefficients that we have been just quoted correspond to massless quark; they have 
to be accordingly modified for heavy quarks. When this is done and comparison with 
experiment is done we obtain from all four LEP experiments 

asiMz) = 0.123 ± 0.008. (66) 

The evolution of a{e'^e~ — > hadrons) as a function of the energy is visualized in fig. 16 



Fig. 16.- a{e^e — > hadrons). From^-*^) 

Some comments are in order here. The first one is that the convergence of the series 
and, ultimately, our ability to get physics out of it hinges on two points. First of all, has 
to be small enough. A coupling constant as just slightly larger produces a contribution 
for the 0{a^g) term comparable to the one of 0{a1), so the convergence of the series 
deteriorates very quickly at low energies. Secondly, the coefficients are actually scheme 
dependent. For instance, if we repeat the calculation but in the MS scheme we get 

^MS = ^MS + (iog47r - 7^)^^I1^ ~ 7.35 - 0.44Ar^. (67) 
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The convergence of the perturbative series is much worse in the MS scheme. It is not clear 
why it is so good in the MS, because in Field Theory this is often the case only if one 
uses a physically motivated scheme (such a subtraction at some energy scale), while the 
reasons to use the MS scheme are of practical order. Be as it may, this is a welcome fact. 

Rhad is an observable so it must be independent of the renormalization scheme. If 
we work in the MS scheme something else must change so that the net result is still the 
same. The thing that changes is of course the coupling constant itself 

^^^""^ ^^^-(33-2iV,)log(^VA|^)log(/.VA|^)' ^''^ 
which just makes up for the changes (up to log log terms), i.e. 

l^^n'^^r-^^^!llM^^^.„ = r + ^+,f^^^f + .„ (69) 

TT TT TT TT 

This is a general feature of perturbation theory. In this framework we necessarily deal 
with truncated series, so independence of the subtraction point or of the scheme can only 
be checked up to terms of the next order in the expansion. That is {R denotes some 
renormalization prescription) 

oo oo 

J2 Cn{R)as{RT = E Cn{R')as{RT, (70) 

n=0 n=0 

but 

N N 

J2 Cn{R)as{Rr E Cn{R')asiRT; (71) 



n=0 n=0 

the error being of 0{a^~^^). 

A lot of theoretical work has been done on the issue of which is the 'optimal' scheme^^) . 
Some possibilities are 

(i) FAC (Fastest Apparent Convergence): choose the scheme that makes cn = 0, cjv 
being the last computed coefficient. 

(ii) PMS (Principle of Minimal Sensitivity): demand that 

d 

——(observable) — 0, 
ait 

and work in the scheme R that fulfills this equation (which, of course, would be any 
scheme if we knew the observable exactly). 

It is a fact that the quality of the series improves when one uses these methods, but 
unfortunately one is forced in general to use different schemes for different observables. On 
the basis of these analysis it has even been claimed that ctg somehow 'freezes' at ~ 0.3 at 
low energies. However, it is fair to say that the general properties of the perturbative series 
in QCD are so poorly understood that any method that does not directly rely on actually 



20 



computing the neglected first term in the perturbative expansion and making sure that it 
is small is likely to be met with skepticism. It is hard to base derivations of fundamental 
parametres such as cug on 'optimization' techniques. 



9.- Rr 

On energy considerations it is obvious that the r is the only lepton we know heavy enough 
to decay into hadrons. This of course makes it a very interesting object. The inclusive 
decay rate r — > + hadrons can, in principle, be derived from QCD by exactly the same 
techniques as Rhad- This is shown in fig. 17. 



Fig. 17.- Determination of Rr through dispersion relations. 
The counterpart of Rhad is now 



^ _ r(r ^ + hadrons) 



At lowest order in QCD 



Rr = nr^r^rdu)+rir-..rsu) ^ ^ 

Rr can be computed as a power series in ctg. Unlike in Rhad, the typical scale will be very 
low (even zero) due to kinematical reasons (see fig. 17). To tackle this problem we decom- 
pose the W boson self-energy (which here plays the role the photon vacuum polarization 
did in Rhad) into vector and axial parts as well as Cabibbo-allowed and Cabibbo-suppressed 
terms 

n'^^ = \Vud\\Kdv + K^iA) + \Vus?{K:v + K^a)- (74) 

Each Hi^^ can be split into transverse and longitudinal parts 

W = i-g'"'k^ + k^'k'')U^^^ + k^'k''U^°^ . (75) 

For massless fermions 11 = 0. We assume that this is the case for simplicity. For a 
detailed account see^^^ . We can write 



Rr = 1277 P ^(1 - 4)2(1 + ^)Imn(i)(.) 
Jo 

= 67r^ / 4(1- 4)2(1 + ^)nW(.) 



(76) 
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where we have used Cauchy theorem to write the integral we are interested in as a contour 
integral. For large |s| we can compute n*^^''(s) just as we did for the hadronic vacuum 
polarization. One gets basically the same result (up to factors). For instance 

Imn(^)(.) = U\Vua\' + |K.|^)[1 + ^ + r,{^f + rsC^f + ...]• (77) 

ZTT TT TT TT 



We also know that 



as{s) _ as{ml) ^ A Q!.(m^) _s_ ^ ^^^^ 



TT TT 2 TT 

So, finally 

Rr = S{\Vud\' + \Vus\')[l + + (r2 - ^ipi){^^r + ... + n.p.t.] (79) 

n.p.t. stands for non-perturbative contributions (contributions non-expressable as a power 
series in ct^). We did not have to worry too much about them in Rhad because there they 
were characteristically suppressed by a power of the momentum transfer. Here they can 
be important. Fortunately, they have been found to be smalP^^. Anyhow, we are now in 
posession of an alternative way of determining ctg though r decay. The best data comes 
again from LEP. A fit to experimental numbers (including mass corrections, which have 
been neglected throughout) gives 

a^(m^) = 0.341 ±0.035 ^^(Mz) = 0.121 ± 0.004. (80) 

Thanks to the logarithmic scaling this method provides us with the most accurate deter- 
mination of as so far. It is also a nice way of testing next to leading scaling. 

10.- Logs in QCD 

Before continuing our discussion, it is convenient to make a short theoretical digression. 
We have seen how the scaling of ag with the energy is logarithmic. Logs in fact play a 
very crucial role in Quantum Field Theory. Let us then stop and think for a second which 
is the origin of these logarithmic terms. 

Actually there are two types of logarithms in a Field Theory such as QCD 

log^ log^ (81) 

They have very difi^erent origin. is some renormalization or subtraction scale, while can 
be some external momentum squared or a small (mass)^ that we have given by hand to the 
gluon. The first type are associated to ultraviolet divergent integrals (integrals with a bad 
behaviour when the internal momentum is large) . The second type are infrared logs and are 
related to Feynman diagrams with a bad behaviour when one or more external momenta 
vanish. Ultraviolet logs appear in any renormalizable Field Theory after renormalization. 
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On the contrary, infrared logs appear whenever a theory has massless particles in the 
spectrum (such as photons or gluons). 

A given Feyman diagram can give rise to both type of singularities at the same time. 
This is illustrated in fig. 18 



Fig. 18.- Infrared logs. 



There actually two classes of infrared logs caused by massless particles. The so-called 
infrared divergences arise from the presence of a soft massless particle {k^ 0). For 
instance in the process e"'"e~ — > jJi^ iJi~ at the one loop level (fig. 19) we have to compute 
the integral 



Fig. 19.- Example of diagram having an infrared divergence. 



' (82) 



(27r)4 /c2[(pi + ky - m?][{j)2 + kf - m? 
When p\= p^ = the integral behaves for A;'^ ^ as 



/ 



d^k 1 



(83) 



and diverges. This divergence is unphysical so it must be cancelled by something else. The 
Bloch-Nordsieck theorem^^) states that in inclusive enough cross-sections the infrared logs 
cancel. What do we mean by 'inclusive enough' ? A detector will not be able to discern a 
'true' muon from a muon accompanied by a soft enough photon (with k ^ 0). Therefore, 
in addition of the diagram shown in fig. 20 we have to consider diagrams where a soft 
photon is radiated by the muon, square the modulus of the amplitude and integrate over 
the available phase space (which actually depends on the experimental cut). When this is 
done the result is infrared finite. The relevant diagrams are depicted in fig. 20 
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Fig. 20.- Real and virtual photons have to be included for IR safe results. 



The other type of infrared logs are called mass singularities. They occur in theories 
with massless particles because two parallel massless particles have an invariant mass equal 
to zero 

k^ = {ki+k2f = ||(a;i+a;2,0,0,a;i + a;2)|| =0. (84) 
The appeareance of such a mass singularity is illustrated in fig. 21 



Fig. 21.- Diagram with a mass singularity. 



{p-kY p2 ^ P - 2A; V + 2A; V cos 6* ' ^ ' 

the denominator vanishes when we set all particles on shell {p^ = /c^ = 0) and ^ ^ (i.e. k 
is parallel to p) . Even if one of the two particles is massive there is a singularity, provided 
the 3-momenta are parallel. 

The Kinoshita-Lee-Nauenberg theorem^^^ ensures that for inclusive enough cros sec- 
tion the mass singularities also cancel. Both for mass singularities and for infrared diver- 
gences there is a tradc-off between A^, the infrared regulator of a massless particle, and 
the energy and angle resolution of the inclusive cross section Ai?, A0. 

In practice, it is better to regulate the infrared logs using dimensional regularization 
(introducing leads to difficulties with gauge invariance). Real gluon emission diagrams 
are regulated by performing the phase space integration in n dimensions. 

There is in fact a lot of physical insight hidden in the infrared logs. We have seen 
that the contribution from the diagram in fig. 19 is infrared divergent, i.e. infinite. Yet, 
physical arguments tell us that the probability of finding a 'bare' isolated muon should be 
zero. We know this because detectors are unable to tell apart a muon from a muon plus one 
soft photon or indeed from a muon plus any number of soft photons. Infrared divergences 
in QED can be summed up and then one sees that the probability of finding an isolated 
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muon is indeed zero and not infinite as the one loop diagram led us to believe. Whenever a 
Feynman diagram is infrared divergent it means that we have forgotten something relevant. 

Let us consider in QED the interaction of a charged fermion with an external source 
and let us expand in the number of virtual photons n (i.e. in the number of loops). The 
total amplitude will be expressed as 

oo 

M{p,p')^J2^n{p,p') (86) 

n=0 

then a calculation shows that 

Mo = mo, 

Ml = rriQaB + mi, 

M2 = mo — h mian + 1712, 



The quantities m„ are IR-finite, while B is IR-divergent. The series in eq. 86 can be 
summed up 

CXD 

M = exp(Q!S) ^ mn, rUn ^ a", (88) 

n=0 

and B can be obtained just from the lowest order diagram. Introducing an IR cut-off A, 
B ~ — logm^/A^, which indeed shows that when we remove the cut-off the probability of 
finding an isolated charged fermion is zero in QED. The addition of soft photons changes 
that result multiplying the total amplitude by a factor ~ (AE/X)'^. There is a trade 
between the infrared regulator and AE, A9. The latter are, of course, detector- dependent 
quantities. 

Although only partial results exist^^\ it is believed that a similar exponentiation takes 
place in QCD. Due to the confinement subtleties it is unclear whether the suppression factor 
is compensated by radiation of soft gluons. Even if this compensation does actually take 
place that would not disprove confinement, only that confinement would have nothing to 
do with the structure of infrared singularities of the theory. 



11.- Jets and ag 

The previous discussion can be summarized in the following way: due to IR singularities 
one is forced to consider cross sections not of individual particles in the final state, but 
rather of bunches of particles, each 'hard' quark and gluon surrounded by a 'soft' cloud of 
gluons and, perhaps, quarks. We will call these bunches 'jets'. 

The Bloch-Nordsieck and Kinoshita-Lee-Nauenberg theorems guarantee the finiteness 
of the cross-sections. We have to define an energy and angle resolution. For instance, if p is 
the momentum of a primary quark (see fig. 22) we can impose that the energy of each soft 
particle in its jet satisfies < epo and also that arg(p, ki) < 5. We will get singularities of 
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the form ctg fog e fog 5 when e,S ^ 0. The specific details depend on the precise definition 
of the jet. Popular jet algorithms have been discussed in this Meeting^^. 



Fig. 22.- Idealization of a jet. 

The situation is unfortunately even more involved because hadrons and not quarks 
and gluons are detected. The evolution of the quarks and gluons produced at high momen- 
tum transfer is perturbative at first, until the average separation of the particles becomes 
0{Aq^j-,). Then the confining potential (and the string picture) takes over. Eventually 
one is forced for anything other than fully inclusive observables (such as Rhad) to introduce 
fragmentation and hadronization models to compute the observable cross-sections. Any 
observable that is not fully inclusive is described by a convolution of two very different 
types of physics 

Observable = Perturbative Non — perturbative 

The perturbative part is, in principle, calculable in QCD as a power series in ag- It is 
affected by some 'theoretical' uncertainties since most observables have been calculated 
up to next to leading order, and not beyond, and higher order effects can be important. 
In addition there is the issue of the choice of an adecuate renormalization scale, which 
sometimes is far from obvious. 

The non-perturbative part has to be modelled. Its relation to QCD and its parametres 
(such as as) is unclear. A considerable amount of cross-checking and experimental feed- 
back is required. In general, the more inclusive the observable is the smaller the unknowns 
coming in from the non-perturbative part are. At LEP there are a number of observables 
that are widely used and where the dependence on the hadronization model is believed to 
be under control such as Rs, thrust or energy-energy correlations. The obvious question 
is do they give consistent results? The question is partially answered in the positive in fig. 
23, showing very good agreement between determinations of ccg from jet topologies. 

A big improvement is obtained after the inclusion of next to leading corrections. They 
are clearly required; the data would be just inconsistent otherwise. Fig. 24 shows the 
energy evolution of Rs, the ratio between three and two jets which is directly proportional 
to as compared to the QCD two-loop prediction. From this observable alone the precision 
on as is approximately ±0.010, which is, by itself, quite remarkable. Not too long ago, at 
PETRA energies, the determination of ag was painstakingly difficult. Why is it so much 
easier at LEP energies than hitherto? 

To understand this we look in some more detail at the thrust of the jets produced. 
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Thrust is defined as 

T = mJ-^^ (89) 

\Pi\ 

ht is the thrust axis, which is varied to maximize T. T takes vahics between 0.5 (spherical) 
to T = 1 (complete alignment). Figure 25 shows the average thrust for a number of 
experiments. The LEP value is close to 0.94, higher than in any previous experiment; 
events are well aligned with the momentum of the primary quark. As a consequence, it is 
much easier to count jets at LEP than in any previous machine. 



Fig. 23.- Determination of cts(M^) from different observables by OPAL^^\ The figure 
includes leading and next to leading corrections and two different hadronization models. 
Observables are (left to right): thrust, oblateness, F-Major, C-Par, i?3. 



Fig. 24.- i?3 for different experiments. 
At LEP energies we are in an energy range where perturbative QCD calculations seem 
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to exhibit a reasonable convergence. If we consider the average value as (// = 34 GeV) = 
0.158 ±0.020 obtained from e+e" machines, excluding LEP (but including Tristan data), 
with very similar, if not identical, hadronization models, jet algorithms, etc. we see that 
an 30% increase in the coupling constant enlarges the experimental error by more that 
a factor 3. Radiative corrections rapidly become hard to control. The size of the error 
is crucial to be able to extract physics from the measurements of a^. For instance, the 
previously quoted value of as leads to = 440^220 MeV, while from the LEP results 
one gets = 250l8o° MeV, a much more stringent result. Incidentally, LEP is the first 

experiment where the 0{al) perturbative contribution is bigger than the error induced by 
hadronization models. 



Fig. 25.- Narrowing of jets with increasing Ec 
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12.- D.I.S.: Free Parton Model 

A brilliant confirmation of the existence of nearly free constituents inside the nucleon was 
provided more than twenty years ago by a series of experiments carried out at SLAC^^^. 
Then it became possible to scatter electrons off nucleons in fixed target experiments with 
a typical momentum transfer ~ 1 — 10 (GeV)^, a kinematical range unexplored until that 
time. The kinematics of Deep Inelastic Scattering processes is shown in fig. 3. 

The virtual intermediate boson is far off its mass-shell and scatters off a quark or 



gluon in a time of 0(1/ y—q^)- Typically quarks and gluons are themselves off-shell by an 
amount of 0{Kqcd)- After the scattering the outgoing particles recombine into hadrons 
in a time of 0{1/ Kqcd)- Thus Deep Inelastic is a two-step process 

(i) Short distance scattering occurs with a large momentum transfer. Well described by 
perturbation theory. 

(ii) Outgoing particles recombine. Not calculable in perturbation theory. However, step 
(ii) can be side-stepped all together for fully inclusive rates. Then perturbation theory 
is adecuate to describe many features of DIS. 

If we place ourselves in the center of mass of the hadron and virtual intermediate 
boson both particles move very fast towards each other. Whatever components the hadron 
contains they will all have moments parallel to P'^, up to transversal motion of 0{Kqcd)- 
Let us write 

= ("90) 
The squared CM energy of the lepton and proton constituent will be 

s={xP+kf ^2xPk^xs. (91) 

We neglect masses (as well as the fact that constituents are off-shell by 0{Kqcd)) The 
final momentum of the constituent is xP -\- q. Therefore 

Oc^(xP + qf ~2xPq + q^, (92) 

so a; = —q^/2Pq. If u is the energy transfer in the LAB system, we can also write 

-q^ 

X = — ^ (93) 
zumN 

niN being the nucleon mass. It is convenient to introduce 

In the lab frame y = v/E and < y < 1. y is thus the relative energy loss of the colliding 
lepton. 

Let us for the time being ignore altogether QCD interactions and let us assume that 
constituents of the nucleons (which we will call partons) are free. DIS will then be described 
by an incoherent sum over elementary processes. The partonic differential cross sections 
in the LAB frame will be 
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• vq, i^qf-scattering 



eq'-scattering 



dcTi^ TTmE 2 , 2.-, n2i /o^^ 

d(^-i7 ,9"^. 2 TrmE r 2 , 2n x2i /n«\ 



^ = Q^^^[l + (l-2/)^]. (97) 



The neutral current sector is dominated by 7 interchange below = M|, so we have not 
bothered to include Z exchange. In eq. 97 Q is the quark electric charge (in units of e) 
and m is the target mass. Since = xP^, we just take m = xtun- Then, for instance, 

^=Q' '""'T"-^ [i+(i-,n. (98) 

GjX ay c[ 

Let u{x)dx, d{x)dx, ... be the number of u, d, ... quarks with momentum fraction be- 
tween X and a; + cZa; in a nucleon. Then xu{x), xd{x), ... will be the fraction of the nucleon 
momentum carried by u,d, ... quarks. We, of course, identify quarks with partons and, 
since we assume that they are free, proceed to sum incoherently over the different scatter- 
ing possibilities. For instance in ep — > eX 

- ^ ' -{uix) +u{x)) + -{d{x) + d{x)) + -{s{x) + s{x))]. (99) 



dxdy s xy'^ 9 9 9 

(We neglect here the possible contribution from the sea of heavy quarks in the nucleon.) 
Other DIS processes weigh differently quarks and antiquarks. For instance, in vp — > fiX if 
-g^ < we have 

d^a ^^^^^2^^^^)^^2^^^^)^-^^^^^_^^2]^ ^^QO) 



dxdy TT 

with Cc = cos 6*0, Sc — sin 6*0, the cosinus and sinus of the Cabibbo angle, respectively. 

The parton distribution functions q{x) are quantities which at present, generally 
speaking, cannot be derived from QCD. The other way round is probably more interest- 
ing: we can learn a lot about the non-perturbativc regime of QCD from DIS experiments 
through these parton distribution functions (PDF). Probably the first thing that one learns 
from them is that gluons are very important. From the SL AC-MIT data^^-* 

Q = U + D + S= f dxx{u{x) + d{x) + s{x)) - O.AA, (101) 
Jo 

Q = U + D + S^[ dxx{u{x) + d{x)+s{x))c^ 0.07. (102) 
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The total fraction of momentum carried by quarks (and antiquarks) is only about 50% 

The rest is carried by gluons, showing that although the naive quark model works very 
well is just a gross simplification as a model of hadrons, at least at large —q^. 

Another example that the quark model fails to describe some basic features of hadrons 
is provided by the 'spin of the proton' problem^) . /x-scattering on polarized targets shows 
that the fraction of the total spin of the proton that can naively be associated to constituent 
quarks is surprisingly small. We shall not dwelve on this matter further here. 

Nevertheless, there are some obvious sum rules for the parton distribution functions 
which can ultimately be explained in terms of the quark model. For the proton 



^0 



dx{u{x) - u{x)) = 2, (103) 



/ dx{d{x) - d{x)) = 1, (104) 
Jo 

I dx{s{x) - s{x)) = 0. (105) 
Jo 

On QCD grounds we expect that this free parton model description of the hadrons becomes 
more and more accurate when — — > oo, — > oo, while keeping x fixed. This limit is 
known as Bjorken scaling and in the strict —q^ = oo limit everything depends just on x. 

The parton distribution functions are actually a function of the energy. At —q^ = oo 
QCD is a free theory, so we can imagine quarks and gluons sharing in equal terms the 
total momentum, in a sort of QCD version of the equipartition theorem of Statistical 
Mechanics. Therefore we expect that in this limit the relation between the momentum 
carried by quarks and the one carried by gluons should be N^Nf /2{N^ ~ !)• This can be 
rigorously justified within QCD^^\ From the above limiting value we see and at higher 
energies the total momentum carried by constituent or valence quarks diminishes and that 
an equally important role is played by particles from the Dirac sea of the nucleon. 

Let us now try to rederive the previous results in a more theoretical setting. Let us 
consider for instance vp scattering. Then 



d{—q^)dv TTS^ 

where 



L>'^ = ^Tr[r{l - - l^W]kakp (107) 

is the trace over the leptonic external lines, and H^i, is given by 

Y,{P\J^.mX{P')){X{P')\J,(Q)\P) = I d4^e^'^^(P|J^(^)J,(0)|P)=Imn^,(g), (108) 



X 

with 



n^.(?) = / d"^e^^^(P|TJ^(^)J.(0)|J'). (109) 
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We decompose H^i^ as 



H^. = -g^,F^ + + —^e^.p^PPq'^F^ (110) 

(If we assume that we are working with non-polarized targets P and q are the only vectors 
at our disposal.) Fi, F2 and F3 are the nucleon structure functions. Using the kinematical 
relations x = —q^/2vmN and y = 2mNi^/ s we get 

d{—q^)dv — usdxdy, (HI) 
[Fi V + F,{1 -y)- Fsxy{l - |)]. (112) 



dxdy 27r 

Let us now compare with the free parton model. We see that (restoring the up index, to 
make apparent that the structure functions are process dependent) 

F^^ix) = cl{u{x) + d{x)) + sl{s{x) + u{x)), 

F^P{x) = 2xcl{u(x) + d(x)) + 2xsl(s(x) + u(x)), (113) 

F^P{x) = 2cl{u{x) - d{x)) + 2sl{-s{x) + u{x)). 

For other processes the actual expressions may vary but the structure functions are always 
linear combinations of the parton distribution functions, i.e. F2{x) = x'^^6iqi{x), etc. 
Note that in the free parton model 

Fl{x) = F2{x) - ^ 0. (114) 

This is the Callan-Gross relation, which actually is not an exact one; it gets modified when 
the q^ dependence is included, i.e. we depart from the strict —q^ = 00 limit. F^ in some 
sense measures the spin of the target. Let us assume for one second that our target has 
spin zero instead of one half. Then we can write with the help of just one form factor 

{X{P')\J,mP)^F{q'){p + p'), (115) 

Plugging this back into eq. 108 we immediately observe that Fi=0. However, experimen- 
tally the Callan-Gross relation is well satisfied, and, on the other hand, F2 is certainly 
non-zero (see fig. 26 below). 

13.- Scaling Violations 

Fig. 26 shows some recent data from the ZEUS collaboration at HERA^°) . 
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Fig. 26.- F2 as measured by ZEUS. 



It is clear from the data that there is some — —q^ dependence in the structure func- 
tions. In other words, there are violations of Bjorken scaling and actually Fi — Fi(x, Q^). 
The free parton model is not completely correct (no big surprise, of course). Our job is to 
try to understand these violations in the framework of QCD. 

Let us recall that all the strong interaction effects in DIS are contained in the hadronic 
tensor H^^. For —q'^ — > cxo the integral in cq. 108 is dominated by the z'^ ^ behaviour, 
so in order to find the q'^-dependence for large values of —q^ we need to find the short 
distance expansion of the product of currents J^(^)J,y(0). We evaluate this expansion in 
three steps^-*^) 

• We decompose 

J^(^)J40) = {d^d,-g^,d^)OL{z,0) 

+ {g^xdpd, + Qp^d^dx - g^xgpud'' - g^,dxdp)Oi''{z, 0) (116) 
+ ... 

O2 is assumed to be symmetric in its indices. The dots stand for terms which are antisym- 
metric in jjLv, which we will not take into account (they would contribute to the structure 
function F3; we will not consider it here). Current conservation is built in this expression, 
which is otherwise completely general. 

• We expand eq. 116 in a complete basis of local operators using the Operator Product 
Expansion (OPE)^^)^ 

Oi'^iz,^) = 5^4,^^^)^'^^ ...^'^"O^i^.^^jO), (117) 
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04^,0) = E<n(^')^'^^ •••^'^"C»iMi...Mn(0)- (118) 

There may be more than one operator with a given set of quantum numbers at a given 
order in the expansions eqs. 117-118 and we have included an index i. The functions 
are called Wilson coefficients. Both sides of eqs. 117-118 must agree when inserted in 
any expectation value or Green function. On dimensional grounds, if the dimension of the 
current is do and the one of is then, as z'^ —>■ 0, the Wilson coefficients 

behave has 

~ (^2)-do + K4(n)-n)+2_ (^^g) 

The combination (io(n) — n is called the twist of the operator. For —q^ oo only ^ 
matters, hence we need to retain the operators of lowest twist, as they provide the most 
singular (hence dominant) behaviour. The contribution of higher twist operators is down 
by inverse powers of = —q^. It is convenient to define the Fourier transform of the 
Wilson coefficients 

on+l r 

^n(Q')-(g2y^^''---^"' = J d^ze'^^z^^...z^-e^{z^). (120) 

We have neglected the tensorial structures that would give less singular contributions in 
the — > limit. Naively the c^((5^) are dimcnsionlcss quantities, pure numbers. If fact 
this is not so. The r.h.s. of eqs. 117-118 is not well defined because when we insert the local 
operators in a Green function we will get additional divergences. We have to renormalize 
them and this changes the Wilson coefficients. As usual, via the counterterms, which in 
perturbation theory are a power series in ctg, a logarithmic dependence in a renormalization 
scale ji appears. On dimensional grounds, the dependence must be through Q^///^. Here 
are the scaling violations we were after. 

• Finally, we need to know the expectation values of the local operators in the 

proton state 

{P\Oi^,...^S^)\P) = A\,nP^.. ■ ■ ■ Pl^r.. (121) 

{P\Olt...,SmP) = Al^^P^P^P,^...P^^. (122) 

This is certainly not the most general decomposition of these expectation values, but other 
tensorial structures will eventually prove to be of higher twist. When we collect all the 
pieces of the calculation and put everything together, we finally get 

FL{n,Q') = f\xx--'FL{x,Q^) = J^^nil^^KniQ^l^''), (123) 
Jo i 

F2{n,Q') ^ f dxx^-^F2{x,Q^) = ^ ^^,n(/^')c^,n(Q', /^')- (124) 
Jo i 

Knowing the Wilson coefficients and the expectation values of the matrix elements we 
can, in the large limit, compute the moments of the structure functions and hence the 
structure functions themselves. The scaling violations come about exclusively through the 
logarithmic dependence on Q'^ j y? of the Wilson coefficients. 
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The product An{iJ>^)cn{Q^ , fJ'^) is //-independent because, as we have just seen, is an 
observable. As befits a renormaUzable theory c„(Q^, satisfies a rcnormahzation-group 
type equation (we assume that there is only one operator of a given dimension and quantum 
numbers, so there is no mixing — we suppress the i superindex for simplicity) 

ifi— + oispQ^^ - loMQ^ i"') = 0. (125) 

The quantity 7o„ is the anomalous dimension of the operator. It is found by determin- 
ing the combination of renormalization constants that makes it finite. Eq. 125 can be 
integrated out 

Cn{Q^^i') = Cn(/.^/.^) ( ^^]^ . (126) 



Leading to following scaling behaviour for the moments of the structure functions 

F,(n,g2)=F^(^,^2)(^^^j ^ ^127) 

which is our final expression. Experiments agree on the whole very nicely with the scaling 
violations predicted by QCD. Taking into accoimt all the subtle points of Quantum Field 
Theory that have gone into the analysis, this provides a beautiful empirical check of the 
theoretical framework 



Fig. 27.- The evolution predicted by QCD confronts experiment. From"^"^^. 



14.- Altarelli-Parisi Equations 

Let us now consider an alternative, and perhaps more appealing, way of deriving the 
evolution equation eq. 127. Let us begin by rewriting it in a differential form. Introducing 
the variable 

1 

t=n^^Sj^, (128) 



35 



we have, for instance, 

dF2{n) lo^asit) 



F2{n). (129) 



dt 47r 

Let us now introduce the Mellin transform of the anomalous dimensions 70^ 



(130) 



Using the decomposition 



F2(a;,t) = a;^%/(a;,t) (131) 
/ 



we arrive at 



dqf{x,t) f'djL i^y^t)PC-)^^cif®P. (132) 

dt Jx y y TT 

These are the AltareUi-Parisi equations^^^ They summarize the rate of change of the 
parton distribution functions with t. 

Let us now try to get a physical picture. We start by looking at the QCD diagrams 
that can contribute to this process. At leading order in we have to consider the con- 
tribution from diagrams (a)-(e) in fig. 31. Since we already know that scaling violations 
are logarithmic we have to look for possible sources of logs. Calculations are simplest in 
an axial gauge. A careful evaluation shows that ultraviolet divergences are absent in this 
case, so they cannot provide us with the logs. Infrared divergences also cancel when all 
diagrams are taken into account (although each one of them is infrared divergent). There 
is only one mass singularity, originating from the real gluon emission diagram (d) that 
survives. So this is our only source of logs. 

We have then to allow that the parton actually carries a different momentum fraction 
y and radiates a coUinear gluon. Since F2, up to some kinematical factors, is just a cross 
section describable in terms of partons we can write 

1 

-F2{x,t) = ^Sf dyqf{y)w2{x,y,t). (133) 

t Jx 



In the free case (no coUinear gluon emitted) we simply have 



y ^ 
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Fig. 31.- QCD diagrams contributing to DIS. 



The mass singularity introduces some t dependence in W2- Using the expression of F2 in 
terms of q{x), 

= 1^ lAy)"-^^^ (135) 

We can iterate the process summing up an infinite ladder of diagrams (fig. 32). This is 
accomplished by replacing qf{y) by qf{y,t); the virginal quark is replaced by a quark that 
has already exchanged many coUinear gluons. 



Fig. 32.- 'Handbag' and ladder diagrams. 

We have been considering F2, but the same procedure can be repeated for any struc- 
ture function. As a simplifying hypothesis we have neglected mixing with other operators. 
In fact, the evolution equation is a {2Nf + 1) x {2Nf + 1) matrix. For flavour singlet 
operators life is more complicated; there is mixing with gluon operators and one must also 
consider gluon parton distribution functions as well 

Q. = — — / —my,t)Pq^g{-)+9{y,t)Pg^g{-)] (136) 

^-^ = '^ f %{y.t)P,^^^^^ (137) 
~ Jx i/ y y 

This is schematically depicted in fig. 33. The detailed form of the Altarelli-Parisi kernels 
at leading order can be found in . Next to leading expressions for them are also available 
in the literature. 
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The analysis of Deep Inelastic Scattering based either on the OPE or on the Altarelli- 

Parisis equations have been amongst the most clear tests of perturbative QCD and the 
best way of determining ag- However, at present the value of Aqcd extracted from LEP 
physics is quite competitive, if not better. 



Fig. 33.- Altarelli-Parisi kernels. 



15.- Parton Distribution Functions 

Although there are some exceptions (notably that of the pion form factor) we do not know 
in general how to compute the parton distribution functions, even for —q^ — > oo. Only 
their evolution can be reliably computed either through the Operator Product Expansion 
of the use of the Altarelli-Parisi equations and this for large enough values of —q'^. 

An interesting issue is the behaviour of the parton distribution functions at the end- 
points X — and x — 1. The large n behaviour of the moments (eq. 123-124) probes the 
X ^ 1 region. Since it is natural to expect that at the kinematical boundaries the parton 
distribution functions vanish, one can make the following ansatz for a; — > 1 

q{x, ^ yl(Q2)(l _ ^yia.iQ'))-! (^3g) 

Demanding that eq. 138 fulfills the q^ singlet evolution equation (136) and (137) leads to 

Likewise, for the gluons we have 

~^0r(2 + .(„.(Q^))) log(l-x) <^*i' 

The constants Aq, A'q and i^q are not calculable on perturbative QCD and depend on the 
specific operator, do is universal. 

When X 1 the gluon distribution functions approach zero more rapidly than the 
quark ones. For large values of x the quark contents of nucleons is the relevant one. 
For small values of x the opposite behaviour takes place, the gluon distribution function 
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eventually becomes dominant. At LHC the cross-section will be greatly dominated by low- 
X physics and the important process there will be gluon-gluon scattering. At the Tevatron 
the quark contents of protons and antiprotons is still dominant. 

To get a semi- quantitative (but not totally satisfactory) picture of the low-x behaviour 
of parton distribution functions, one inverts eqs. 123-124 via an inverse Mellin transform. 
Then 

F2ix,t) = ^-j dnx-^''-'-^ F2M (142) 
The evolution of F2(n, t) is known. We can write 

F2M = exp[/(n,t)]F2(n,to) (143) 

Then we proceed to evaluate the r.h.s of eq. 142 by the saddle point method. Of course 
this requires that log(l/a;) is large and this is the reason why this procedure gives only 
the small x behaviour. The solution is parametrized in terms of ^2(^^-0)^0)) "-o being the 
solution of the saddle point equation for n. Working things out one finds that indeed the 
gluon parton distribution function is dominant for low x behaving as 



g{x)r.-expJc{Q^)\og- (144) 

X \ X 

where C{Q'^) is calculable. Unfortunately, this answer is not totally satisfactory beacuse 
something must stop the growth in g{x) for low x, or else one runs into unitarity problems 
sooner or later, and thus eq. 144 it is not credible all the way to a; = 0. Technically 
speaking, there must be corrections that destabilize the saddle point solution. Physically, 
the uncontrolled growth of the gluon distribution is an infrared unstability. The density of 
soft gluons is too large. Shadowing and non-linear evolution equations are the buzzwords 
here^^) . 

Low X means —q^ large, but fixed, and — > 00. This is the Regge limit and we 
understand very little of QCD is this regime. There are several parametrizations*^^) of the 
quark and gluon distribution functions that fit the experimental results rather well (~ 1%) 
in the range 1 > x > 10~^. These parametrizations are obtained in the following way. One 
fits the experimental data at some relatively low value of using, for instance, a form 
for the parton distribution function inspired in the quark model, say q{x) = Ax°'{l — x)^. 
Then one evolves q{x) using the Altarelli-Parisi equations and performs a global fit. The 
region below 10~^ had not been explored experimentally until very recently; a first look 
at these low-x values has been provided by the commissioning of HERA. 

HERA is a machine ideally suited for an in-depth analysis of structure functions. 
Fig. 34 shows the physics reach of HERA. It should be possible to arrive at very low 
values of x (up to a; ~ 10~^). There are already very interesting results on the region 
10""^ < X < 10~^, which actually show that most of the parametrizations of the structure 
functions perform very poorly when it comes to reproducing the data at low-x. Typically 
they predict an increase as x ^ which is lower than what is actually seen (fig. 35). The 
behaviour ^2(0;) ~ x~^, with A~l/2asa;— >0, which is predicted from the Kuraev- 
Fadin-Lipatov evolution equation^^) seems to stand the comparison with HERA results 
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best. However, this behaviour is still incompatible with unitarity and cannot hold all the 
way to a; = either. Obviously there are many open questions in this area. 



Fig. 34.- The physics reach of HERA. From^'''). 



Fig. 35.- Low X structure functions at HERA^^^. 

Complementary information on the parton distribution functions is provided by the 
use of sum rules. If one of the operators that appear in the OPE corresponds to a conserved 
current such as, for instance, the non-singlet operator 

u-ff^u - djf^d. (145) 

The corresponding anomalous dimension vanishes and its expectation value corresponds 
to a conserved charge measured in the nucleon state^^^ . Then we know 

/<i..-F<-.W = ^..> (146) 
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for non-singlet operators and a linear combination of 



/^F<^)(.) (147) 

for the singlet operators. In this case the conserved current is the energy-momentum 
tensor. As an example let us consider ep scattering. The non-singlet conserved current is 
V'7mQV' {Q being the charge matrix) and the associated conserved charge is the nucleon 
electric charge. Then 

J dxx-^FlP{x) = ^(1 + 0{as)). (148) 

Other sum rules are 

/...-(^7W-^rW)=0. (149) 
due to Adler, or the Gross-Llewellyn-Smith sum rule (/: isoscalar target) 

J dxx-'^F^^{x) = 3(1 + 0{as)). (150) 



16.- Large Rapidity Gaps 

It would not be appropriate to conclude this review without saying a few words about 
some very interesting events that have been seen at HERA^^^. 

We discussed in section 7 the phenomenon of color coherence. In DIS we expect 
that hadron production should take place predominantly in the color string joining the 
struck parton and the nucleon remains. In practice this means that most hadron energy 
deposition should occur in an angular region relatively close to the proton beam direction. 
Described in terms of the pseudorapidity rj — — log tan 6*72 (6* being the polar angle mea- 
sured following the incoming proton) one should expect hadron energy deposition to be at 
relatively large values of rj (rjmax = 4.3). 

While this is the case for the bulk of the event sample (fig. 36), a small but sizeable 
fraction remains all the way to ~ —2, with little hadronic activity in between. These 
'large rapidity gap' events cannot be accounted for in the parton model wc have discussed. 

It has been suggested that the large rapidity gap events can be explained in terms of 
the pomeron^-*^) , which is believed to dominate elactic and driffactive production in hadron- 
hadron interactions. Although the pomeron itself falls completely outside perturbative 
QCD (presumably you have to sum infinite families of diagrams to reproduce pomeron- 
like behaviour), the concept of pomeron structure functions has been put forward and 
studied. The idea is that the pomeron could exhibit a partonic structure that could be 
probed in hard diffractive processes. The events observed ar HERA appear to be of this 
kind. 



41 



Fig. 36.- Rapidity distribution of events 



17.- Summary and Outlook 

QCD is now twenty years old. Obviously we are a long way from the days when it was 
hotly debated whether a sensible Field Theory could possibly be asymptotically free. We 
know today that QCD makes perfect sense, in fact QCD is now the paradigm of a Quantum 
Field Theory and the primary suspect is QED, which is believed not to be a consistent 
theory all the way to zero distances. This is indeed a paradoxical situation. 

Quantum Chromodynamics is a strongly coupled theory. Ultimately this can be traced 
to the fact that Aqqd is similar in magnitude to the light quark masses. For very heavy 
quarks perturbation theory gives reasonable estimates, but fails for light quarks where 
confinement plays a crucial role. 

Then it would have seemed rather hopeless twenty years ago to expect that we would 
be able to measure as with a mere 5% error. Confinement makes most hadronic observables 
to depend on as in an unknown form. Non-perturbative studies such as lattice QCD can, 
in principle, unveil the as dependence, but we are still a long way from reaching the 
above precision, even though these studies provide overwhelming evidence that QCD does 
describe the hadronic world. 

Although many valuable theoretical developments exist by now, it is through experi- 
ment that we have ultimately learnt most of what we know about QCD and as- At LEP we 
are in the privileged situation of having an experiment with high statistics, a clean set-up 
and typical processes with a large momentum transfer, leading to a reasonable convergence 
of perturbative expansions for inclusive processes. Thanks to this triple conjunction we 
have been able to reach the level of accuracy of a few per cent. 

Yet, with practically everything computed at the next-to-leading order and with no 
other experiment of similar quality in sight, it is hard to imagine a sizeable reduction in 
the error of as in the near future. 

HERA will probably capture most of the interest of workers in the field in years to 
come. There we have the possibility of exploring a new kinematic range in Deep Inelastic 
Scattering. Sensitivity to the low-x range in structure functions will allow us to explore 
the region where parton wave functions begin to overlap and perturbation theory fails to 
describe the physics even when ^ oo. 

A spectacular evidence for the insufficiency of perturbation theory even for processes 
with large is to be seen in the detection of large rapidity gap events, interpreted in 
terms of collective modes of dual theory. This is not really new; Regge theory has been 
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long required to understand the elastic and low pt behaviour of hadron interactions. What 

is new and interesting is that at HERA we will probably be able to interpolate smoothly 
between the Regge limit of QCD and ordinary perturbation theory. Many surprises are 
probably in store for us. 
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